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Possible direct measurement of the expansion rate of the universe 
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A new method is proposed for directly measuring the expansion rate of the universe through very 
precise measurement of the flux of extremely stable sources. The method is based on the definition 
of the luminosity distance and its change along the time due to cosmic expansion. It is argued that 
galaxies may be chosen as the targets of the observation to perform such measurements. Once the 
required precision of the flux is achieved, profiting from the abundance of galaxies in the universe, 
the method could be quite promising. 
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I. INTRODUCTION 

Measuring the expansion rate of the universe along the 
redshift has been one of the most important scientific ob- 
jectives in cosmology since the discovery of the expansion. 
It is usually pursued by measuring the distances at dif- 
ferent redshifts, which, in turn, could be done with the 
data of the standard candles like the type la supernovae 
(SNe la) Q, Gamma-ray bursts (GRBs) 0, Q etc or 
the standard rulers from cosmic microwave background 
(CMB) and baryon acoustic oscillations (BAO) 5]. 
In fact, it is the measurement of the luminosity distances 
of SNe la [1, 0] that leads to the discovery of the unex- 
pected accelerating universe, which is attributed to the 
mysterious dark energy. The further study of the na- 
ture of the dark energy requires more precise expansion 
history of the universe. The ongoing joint dark energy 
mission (JDEM) |^ has aimed at such a precise mea- 
surement. Currently, the expansion rate of the universe is 
still mainly measured through the distance measurement. 
The distances depend on the expansion rate through an 
integration, so the extraction of the expansion rate from 
the distances involves differentials, which significantly af- 
fects the precision of the measurement of the expansion, 
not to mention that the derivation of the dark energy 
equation of state (EOS) from the expansion rate involves 
differentials again. 

On the other hand, despite of the difficulties, some 
proposals have been presented for directly measuring the 
expansion rate of the universe, for example, through the 
measurement of radial BAO Q , the relative ages of pas- 
sively evolving galaxies (loj . the temperature and po- 
larization anisotropies of the CMB [11[, or the redshift 
drift [12, 13 (the socalled Sandage-Loeb test). In this 
paper, we propose another method to directly measure 
the expansion rate by measuring the changes of the flux 
of extremely stable sources. 
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II. METHODOLOGY AND DISCUSSION 

Consider a source rest at comoving distance r, with a 
redshift of z, its luminosity distance to us is given by 



a[to)r{l + z), 



(1) 



where a is the scale factor as a function of time and to 
denotes the time of today. We assume the signal we ob- 
served at the time of is emitted by the source at tom- 
Due to the expansion of the universe, if we observe the 
source again after a time interval of A^o, i-e. at the time of 
io+A^o (the corresponding signal is emitted by the source 
at the time of icm + Atcm with Aicm = Atp/ (1 + -z)), we 
will find its luminosity distance changed with a value (we 
only take into account the first order terms in this paper) 
of 



Adi = /\a{to)r{l + z) + a{to)rAz 



Aa(to) 
a{to) 



Az 



l + z' 



(2) 



where Aa(to) and Az are the changes of a(io) and z in the 
time interval Ato due to the expansion of the universe. 
We can rewrite Eq. ([2]) into 



Adr 



Aa(io) 



a{to) l + z 
For the expansion of the universe, we have 
Aa(to) a(io) 



(3) 



a(to) 



a{to) 
HoAto 



Atn 



(4) 



and since 1 + z = a{t^)/a{t. 
Az = 



a(^o) gfa) d(^cm) A^ 
^lo rAtn 



a(iom) a(iem) a(tem) 

4H^Ato-(l + z)i/(.) 
[(1 + z)Ho - H{z)] Ato, 



Ato 

l + z 



(5) 
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where Hq and H{z) are the Hubble parameter of today 
and that at the redshift of z correspondingly. Eq. ([SJ is 
in fact the core of the Sandage-Loeb test O [3, i-e., if 
we manage to measure the redshift drift Az, we obtain 
the Hubble parameter. Substitute Eq. (jj]) and Eq. ^ 
into Eq. ([3]), we have 



Adz 



2Ho 



H{z) 



l + z 



Ato. 



(6) 



With this equation at hand, one may naturally think 
that, similarly as the Sandage-Loeb test, if the luminos- 
ity distances can be measured to a very high precision 
such that we could distinguish the small changes in the 
luminosity distances for a reasonable time duration, we 
could also immediately derive the corresponding Hubble 
parameters. Unfortunately, the distance measurement it- 
self is a difficult task in astronomy, especially for cosmic 
distances. One of the most precise ways of measuring 
cosmic distances is through the observation of gravita- 
tional waves, for which a relative precision of about 10"'^ 
is expected for the luminosity distance |14|. From Fig- 



ure [1] we can see that, even for such precise measure- 
ments, about 10* years will be needed before we could 
tell the changes in the luminosity distances, which is ob- 
viously not feasible. 




FIG. 1. AdL/dL per year versus redshift z for the flat 
ACDM cosmological model with Om,o = 0.27 and Ho = 
70.5 kms"^ Mpc-^ 

Recall that the luminosity distance is defined through 

L 



F = 



(7) 



where L is the luminosity of the source and F is the 
observed flux, we have 



AF _ AL ^l^dL 



(8) 



where AF and AL are the changes in F and L during 
the time interval Ato- So, if we point our telescope to 



some extremely stable source such that AL/L can be ig- 
nored in Eq. and manage to measure the flux F to 
a very high precision, we could measure l^d^jdh indi- 
rectly through AF/F. Thus, Eq. (g]) could stiU be used 
to measure the expansion rate of the universe. Then the 
problem becomes whether there exist such extremely sta- 
ble sources and whether we could manage to measure the 
flux F to the precision needed. 

For the requirement of extremely stable sources, we 
may consider objects which include lots of similar 
sources, so that we could statistically reduce AL/F to 
a very low level. For example, a galaxy includes lots of 
stars. During the time interval A^o, the luminosity of a 
star in the galaxy may increase or decrease. Let Li be 
the luminosity of the ith star in the galaxy, we may view 
l^Li/Li as a random variable and, for a simple estima- 
tion, assume it follow the normal distribution A/'(0, u^), 
then, for the galaxy. 



AF 
17 



AA(0, 



(9) 



and since J2i ^1 /iJ2i ^i)'^ has the order of 1/N, where 
N is the total number of stars in the galaxy, the stan- 
dard deviation of AL/L for the galaxy is thus reduced 
by a factor of about 1/VN compared to that for a star. 
While considering the change in the luminosity, we first 
exclude galaxies involving violent astrophysical processes 
like stellar explosions, which usually could be easily iden- 
tified. If we simply model the luminosity evolution of a 
steadily burning star to be of that it linearly increases 
to its maximum value and then linearly decreases, then 
AL/L for the star should have the order of the inverse 
of its lifetime. The lifetime of a star ranges from only 
a few million years (for the most massive) to trillions of 
years (for the least massive) . Here, we conservatively set 
the AL/L for stars to be of the order of 10~^ per year 
which corresponds to stars with smallest lifetime. Typ- 
ical galaxies consist of from 10'' to 10^^ stars. If we set 
N — 10^°, the standard deviation of AL/L for the galaxy 
would be of the order of 10^^^ per year. So the AL/L 
can hopefully be reduced to below the expected values 
of Ad^/dL- In practice, some selection criteria on the 
galaxies may be needed. The details on the selection are 
out of the scope of this short letter. But it is worth men- 
tioning that one may concern the impact from the evolu- 
tion of the observed galaxy beyond its stability. On this 
issue, first we should, of course, select passively evolving 
galaxies as our targets of observation basing our knowl- 
edge of galaxies themselves. Second, the AF/F caused 
by the cosmic expansion has a flat-line spectrum, while 
those caused by other astrophysical processes usually do 
not has such a character. With this, we could further 
exclude those galaxies whose evolution (including effects 
from star-formation, dust, etc.) dominates over the cos- 
mic expansion on AF/F. 

For the measurement of the flux, cryogenic detectors 
usually can achieve very high precisions. See, for exam- 
ple |15l4l7l |. for cryogenic detectors. We can also improve 
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the precision by time integration, and/or fitting the flux 
observations along the time to a linear relation. Cross- 
correlation between observations of different bands may 
also help. The requirement on sensitivity of the instru- 
ment for the observation, when taking all the considera- 
tions into account, need further investigations. 

Anyway, if we do manage to measure the flux to a very 
high precision, profiting from the abundance of galaxies 
in the universe, this method may be quite promising for 
measuring the expansion rate of the universe. Observa- 
tion of galaxies at the same redshift may be averaged 
to increase the accuracy of the measurement of the ex- 
pansion rate, or hopefully, if the accuracy achieved with 
single galaxies is good enough, we could even map the 
three-dimensional cosmic expansion, i.e. the cosmic ex- 
pansion of the different directions in the sky along the 
redshift. Also note that Adi/di has the same order 
from the redshift of zero to redshifts z > 4 (see Fig- 
ure [I}, where the dark energy plays its role in the cosmic 
expansion, so the method may track the whole dynamics 
of the dark energy. 

Compared to the Sandage-Loeb test, though the 
derivation of our method is very similar to that of 
Sandage-Loeb test, the requirement on the instrument, 
the observational target, and the expected outcome are 
quite different. The Sandage-Loeb test requires precise 
measurement of the redshift, while our method requires 
precise measurement of the flux. The targets for the 
Sandage-Loeb test are quasars, while in our method, 
the targets are galaxies, which are much more abundant 
in our universe. Provided the measurement from sin- 
gle galaxies could achieve an accuracy good enough, our 
method would give more detailed information about the 
cosmic expansion across the sky. The redshift coverages 
of the two methods are also quite different. As shown 
in [l^, the redshift coverage of the Sandage-Loeb test is 
roughly between 2 and 5, while our method could poten- 
tially cover the range from a redshift close to zero (see 
discussions in the next section) to redshifts z > A. All 
these make our method a novel one from the Sandage- 
Loeb test despite the similar derivation. 

In the above analyses, we have only considered the 
ideal condition, for example, we have assumed the source 
rest in the comoving reference frame and the photons 
propagate freely from source to the observer. This is 
of course not the case of that in the real universe. For 
the very high precision required by the method, impacts 
from the peculiar velocity of the source and from the 
gravitational lensing should be investigated seriously. In 
the next section, we showed that the impacts from the 
peculiar velocity can be safely ignored in a wide redshift 
range. For the gravitational lensing, since we have so 
many galaxies in our universe, we can expect that its 
impacts could, at least, be eliminated by an average over 
the galaxies at the same redshift. We leave it to the future 
studies whether the gravitational lensing will affect the 
measurement of single galaxies in our method. 



III. PECULIAR VELOCITY 

In this section, we studied the impacts from the pecu- 
liar velocity of the source on our method. 

First, consider a isotropic point source and an observer 
nearby rest in a local inertial reference frame. Assume 
the source have a velocity of v with respect to the ob- 
server. Imagine a spherical coordinate with the source 
as the origin and the direction of v as the zenith direc- 
tion. Say, the observer have an inclination angle a in 
the coordinate. The momenta of the photons received 
by the observer would have the same inclination angle a. 
But, to the source, as a result of the aberration of light 
caused by the motion of the source, the momenta have a 
different inclination angle a, which relates a through 



(10) 



where /3 — v/c. So, the corresponding solid angles around 
above mentioned momenta have different magnitudes to 
the observer and to the source. The ratio between them 
is 




sin a da 



J{a,(3). 



(11) 



sin a da 

In addition to this, there is also the Doppler effect 

1 + {1- (3 cos a)7, (12) 



where Zd is the Doppler redshift and 7 = 1/^1 — f3^. 
Taking into account both the aberration of light and the 
Doppler effect, compared to the case of a rest source, the 
flux observed by the observer is increased by a factor of 



1 



1 



1 



(l + zrf)2 J(a,/3) [(l-/3cosa)7]2J(a,/3)' 



(13) 



Next, consider the cosmological case of a isotropic 
point source at comoving distance r, with a luminosity 
of L, a peculiar velocity of v, and a cosmological redshift 
of Zc, and the direction of v is at an angle of 9 relative to 
the line of sight (from the observer to the source). The 
photons emitted by the source first experience an aber- 
ration and the Doppler effect due to the velocity of the 
source, then the cosmological redshift before they reach 
us. As a result, the observed flux becomes 



1 



1 



1 



1 



inalr^ (1 + z,)^ (1 + Zd)^ J{tt - 9, (3) 

1 1 1 L 

Analr^ (1 + Zcf [(1 + /3 cos 9)^]"^ J{7t -9,I3) ' 

(14) 

where oq = a(io)- So, when taking into account the 
peculiar velocity of the source, the corrected luminosity 
distance is given by 

dL,p = aor{l + Zc){l + Zd)\/ J{n - 9,13) 



a^r {\ + zc){l+f3 cos 9)-t^ J {tt - 0,/3). (15) 
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It is easy to check that, under the condition of , 
Eq. (dni) and Eq. ^ reduce to 



J(a,/3) 



and Eq. ([T5|) reduces to 



a ~ a + /Jsina, 
1 



(1 + /?C0SQ;)2 



«1, 

(16) 
(17) 



dL,p = aor{l + Zc){l + Zd)/(1 - [3 cos 9) 

= aor(l + Zc)(l + /3 cos 61)7/(1 - /3 cos 6*). (18) 

Further ignore the change in the velocity of the source 
during the observation, we have 



Ad 



Ar 
r 



1 + Zc 



(19) 



We plot - ^) along the redshift in Fig- 

ure [21 from which we can see that, from very low redshift 
to very high redshift, the impact of the peculiar velocity 
is very small and can be safely ignored. 




Comparing the right hand side of this equation with that 
of Eq. (Ul, we can see that the first term is unchanged. 
The second term is newly introduced corresponding to 
the distance the source moved through during the obser- 
vation. While the last term, at first glance, has the same 
form as the last term of Eq. ([3]) , it actually includes two 
parts, which may be called the time part and the space 
part. The time part arises from that the redshift of a rest 
source changes along the time due to the cosmic expan- 
sion, as was shown in Eq. ([5]). The space part is in fact 
that the source has moved to a different location during 
the observation due to its peculiar velocity and sources 
at different comoving distances have different redshifts. 
It is given by 



Az 



l + Zc 



_ H{zc) apAr 

p ^ c(l + zc) TT^l^' 



(20) 



where the subscript p denotes the redshift drift caused by 
the peculiar velocity, i.e., the space part. So, the impact 
of the peculiar velocity on the change of the luminosity 
distance is 



Ad 



L,p 



AdL 

dL,p dL 
where 



H{z,) 



c(l + zc)2 Vl - kr^ 



Ar 



(21) 



Ar _ (1 + Zcf^/1 - kr'^ a{tc„,)Ar 
r 



dL 

\2. 



(1 + Zc) Vl - fcr-2 



dL 



Vl — kr"^ 
V cos dAt„ 



(1 + Zc)Vl - kr^ 
dL 



V cos 9Ato 



and 



\/ 1 — kr^ 




kfi- 



Hi 



c2 {l + z,Y 



(22) 



(23) 



FIG. 2. - ^) versus redshift z for a peculiar 

velocity oi v = 1000 kms~^ and the flat ACDM cosmological 
model with f2,„,o = 0.27 and Hq = 70.5 kms~^ Mpc~^. 



IV. SUMMARY 

In summary, basing on the definition of the luminosity 
distance and its change along the time due to cosmic ex- 
pansion, we discussed the possibility of directly measur- 
ing the expansion rate of the universe through very pre- 
cise measurement of the fiux. Extremely stable sources 
are needed for the method. We argued galaxies may be 
chosen as the targets of the observation. Since composed 
of many stars, their luminosity is very stable. Those 
involving violent astrophysical processes during the ob- 
servation could be easily identified and excluded. Fur- 
thermore, the relative flux change caused by the cosmic 
expansion has a fiat-line spectrum, while those caused 
by other astrophysical processes usually do not has such 
a character. With this, we could further exclude those 
galaxies whose evolution dominates over the cosmic ex- 
pansion on the flux change. We also showed that the 
peculiar velocity can be safely ignored for a wide redshift 
range in our method. Once the required precision of the 
fiux is achieved, profiting from the abundance of galaxies 
in the universe, the method could be quite promising. 
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